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We show the existence of a statistically significant, robust detection of a gamma-ray source in the 
Milky Way Galactic Center that is consistent with a spatially extended signal using about 4 years of 
Fermi-LAT data. The gamma-ray flux is consistent with annihilation of dark matter particles with 
a thermal annihilation cross-section if the spatial distribution of dark matter particles is similar to 
the predictions of dark matter only simulations. We find statistically significant detections of an 
extended source with gamma-ray spectrum that is consistent with dark matter particle masses of 
approximately 10 GeV to 1 TeV annihilating to bb quarks, and masses approximately 10 GeV to 
30 GeV annihilating to rf leptons. However, a part of the allowed region in this interpretation 
is in conflict with constraints from Fermi observations of the Milky Way satellites. The biggest 
improvement over the fit including just the point sources is obtained for a 30 GeV dark matter 
particle annihilating to bb quarks. The gamma-ray intensity and spectrum are also well fit with 
emission from a millisecond pulsar (MSP) population following a density profile like that of low- 
mass X-ray binaries observed in M31. The greatest goodness-of-fit of the extended emission is with 
spectra consistent with known astrophysical sources like MSPs in globular clusters or cosmic ray 
bremsstrahlung on molecular gas. Therefore, we conclude that the bulk of the emission is likely 
from an unresolved or spatially extended astrophysical source. However, the interesting possibility 
of all or part of the extended emission being from dark matter annihilation cannot be excluded at 
present. 

PACS numbers: 95.35. +d,95.55.Ka,95.85.Pw,97.60.Gb 



I. INTRODUCTION 

The successful launch and operation of the Large Area 
Telescope (LAT) aboard the Fermi Gamma-Ray Space 
Telescope has mapped the gamma-ray sky with unprece- 
dented precision [1 j. One of the principle scientific ob- 
jectives of the Fermi-LAT is to probe the nature of dark 
matter [2], since the canonical weakly interacting mas- 
sive particle candidates' thermal production process in 
the early universe requires significant annihilation in dark 
matter overdensities today if the dominant annihilation 
channel is s-wave [3]. Numerical studies that do not in- 
clude star formation have found that cold dark matter 
particles have a density profile that is strongly centrally- 
peaked [ij [5] . This leads to a galactic gamma-ray lumi- 
nosity from dark matter annihilation that also strongly 
peaks at the Galactic Center (GC) [6 . The largest lu- 
minosity signal arises from the Milky Way Galactic halo 
itself instead of un-associated halo substructure or extra- 
galactic sources [7]. Tempering this optimistic outlook 
for dark matter detection is the fact that the GC also 
harbors a large number of astrophysical sources with a 
high integrated luminosity in gamma-rays. 

Results from observations of the 3° x 3° region about 
the GC by Fermi-LAT have placed competitive con- 
straints on annihilating dark matter (e.g., Ref. [8]). How- 
ever the best, robust constraints on annihilating dark 
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matter come from the much lower-background stacked 
observations toward the dark matter halos associated 
with dwarf galaxies [9, 10 . Observations toward the 
GC by the High Energy Stereoscopic System (HESS) 
telescope have the greatest sensitivity to dark matter 
above a dark matter particle mass « 500 GeV and place 
the strongest constraints on annihilating dark matter 
above that mass [TTJ [12]. This is primarily because 
astrophysical backgrounds are largely reduced at these 
higher energies once the signal from the Galactic Ridge 
is masked [13] . There has also been a set of analyses 
of the Fermi-LAT data toward the GC that find a sig- 
nal consistent in morphology and spectrum with roughly 
10 — 40 GeV mass Weakly Interacting Massive Particles 
(WIMPs) annihilating into r leptons, b quarks or a com- 
bination of both channels [T4HT6] . The spectrum and 
amplitude of the signal was shown to also be consistent 
with a population of millisecond pulsars in the Galactic 
Central stellar cluster [17 and radiation from cosmic ray 
interaction with gas in the GC region [T6 ] H8 ] [T9], 

Pioneering work using EGRET data had previously 
found that emission from the GC could be consistent 
with a WIMP with roughly thermal annihilation cross 
section and ~50 — 500 GeV particle mass, and had also 
forecast that Fermi-LAT would be able to resolve the 
spatial extent of the structure [20 . Preliminary analyses 
by the Fermi-LAT Collaboration did not report evidence 
of an extended source in the GC, though an excess in 
observed counts to model is seen in their results near en- 
ergies of 2 — 5 GeV [21] . In another independent analysis 
of the GC using photons from 1 - 300 GeV, Ref. [22] 
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found that the log-likelihood improved considerably (25) 
with an additional component that had the same spatial 
morphology as that in Ref . [TJi • There has also been con- 
siderable interest in evidence for a line signal associated 
with the GC [23 . The regions used for the line signal 
include a larger area on the sky than what is evaluated 
here, and so we do not discuss that aspect of dark matter 
annihilation signal here. 

In this paper, we present the analysis of 3.8 years of 
data from the Fermi-LAT in the inner 7° x 7° toward 
the Milky Way Galactic Center using the current second 
year Fermi-LAT point source catalog (2FGL), the second- 
year Fermi-LAT diffuse Galactic map, isotropic emission 
model, and new models for any extended emission coming 
from the GC. We find that due to the required fitting of 
the point sources and known extended sources with any 
new sources, there exists a degeneracy between the spec- 
tral properties of point source emitters in the inner < 1°, 
the Galactic diffuse model, and any new extended source 
in the GC. Despite this degeneracy, we find that there is 
a statistically significant, robust detection of an extended 
source not present in the 2FGL or diffuse Galactic map 
that can be consistent with astrophysical or dark mat- 
ter annihilation sources. We discuss both possibilities in 
detail below. 



II. MODELS FOR SPATIAL MAPS OF 
EXTENDED SOURCE 

The excess emission from the GC is centrally-peaked, 
so we only consider centrally-peaked dark matter halo 
models. The dark matter halo models we include in this 
study are the "a/37" profiles fashioned after the Navarro- 
Frenk- White (NFW) profiles [H[25], 



p(r) = 



(r/r s ) 7 (1 + (r/r E 



\a\(/3-7)/a 



(2.1) 



with fixed halo parameters a = 1, f3 = 3, r s = 23.1 kpc, 
and a varied 7 inner profile. The canonical NFW pro- 
file has 7 = 1. Higher resolution simulations show that 
the inner log-slope does not asymptote to a constant 
but rather becomes softer. We also include a fit to 
an "Einasto" profile because higher-resolution numeri- 
cal dark-matter only simulations seem to prefer this fit 
where the log-slope rolls with decreasing radius [26] [27] , 



PEinastoO) = p s exp 
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(2.2) 



with o;e = 0.17 and r s = 20 kpc. This Einasto pro- 
file should be considered as being more representative 
of the DM-only simulations. Substantially more peaky 
profiles like the 7 = 1.2 require other physics, such as 
baryon-induced adiabatic contraction of the halos, e.g., 
see Ref. [28 . However, recent simulations also go the 
other way in that the feedback from supernovae reduce 
the density of dark matter in the center [29] . Note that 



the differences between a 7 = 1.2 and Einasto profile 
(with the parameters as fixed above) are about a factor 
of 2 in the inner 100 pc (or about 0.7 degree) and about 
a factor of 5 in the inner 10 pc. Since the annihilation 
flux goes as density squared, these are substantial differ- 
ences. We also discuss later that the annihilation flux 
maps resulting these peaked density profiles may also be 
the appropriate distribution for an unresolved MSP pop- 
ulation in the inner Galactic bulge region. 

The differential flux for a dark matter candidate with 
cross-section (<jav) toward Galactic coordinates (b,£) is 



d$(M) ((TAv)J(bJ) 1 dN 7 



dE 



Jo 47rm| dE 



(2.3) 



where dN^/dE is the gamma-ray spectrum per annihi- 
lation and m x is the dark matter particle mass. The 
quantity J is the integrated mass density squared along 
line-of-sight, x, 



(2.4) 



J(M) = Jo J dx p 2 (r gal (b,£,x)) 
where distance from the GC is given by 



r gal (b,£,x) = \JR%- ZxRq cos(f) cos(6) + x 2 . (2.5) 

Here, Jo = l/[8.5 kpc(0.3 GeV cm -3 ) 2 ] is a normaliza- 
tion that makes J unitless and cancels in final expres- 
sions for observables. The value for the solar distance is 
taken to be Rq = 8.25 kpc [50] . The density p s in both 
the a/3j and Einasto profiles is a normalizing constant 
degenerate with the local dark matter density, p . We 
adopt a conservative (broad) range of local dark mat- 
ter densities consistent with its most robust estimates, 
Pq = 0.3 ±0.1 GeV cm -3 [501 EI]- However, the spa- 
tial profile on the sky is independent of this uncertainty 
and only is relevant when converting from a flux to the 
particle annihilation rate, as discussed in ^V] 

We also consider flux maps that are proportional to 
projected density profiles as is appropriate, for example, 
when the extended source is the result of the superposi- 
tion of unresolved sources. In this case, 



dE 



dE 



(2.6) 



where F(b,£) is normalized to unity within the region 
of interest and the flux normalization is included in the 
spectrum dN 1 /dE. We consider both spherically sym- 
metric and axisymmetric models such that, 

F(M)oc j dx(r gal , a (M^)r r_1 , (2.7) 

(r g ai,a(M,20) 2 = (Rq ~xcos(b) cos(f)) 2 
+(x cos(b) sin(f)) 2 + (x sin(6) /a) 2 . 

For spherical symmetry (a = 1), the flux map in the 
central parts (of interest here) can be approximated (to 
about 10%) as F(b,£) ex l/(l 2 + b 2 ) T/2 - 
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FIG. 1. Shown in the top row are photon counts in four energy bins that have significant evidence for an extended source 
with a spectrum, morphology, and rate consistent with a 30 GeV mass WIMP annihilating to 66-quarks in the 7° x 7° region 
about the GC. This row shows the 17 2FGL point sources in the ROI as circles. The second row shows the residuals for the 
fit to the region varying all the sources in the 2FGL catalog as well as the amplitudes of Galactic diffuse and isotropic diffuse 
models. The presence of an extended source and oversubraction of the central point sources are visible here. The third row 
shows the best fit model counts for 30 GeV WIMP annihilating to 66-quarks. The fourth row is the residual emission for this 
model without subtracting the extended component. The fifth row contains the residuals when the extended component is also 
subtracted. The maps have been filtered with a Gaussian of width a — 0.3°. 
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III. METHOD 

We use Fermi-LAT data from August 4, 2008 to June 
6, 2012 in the extended source analysis, downloaded from 
the LAT data server at the Fermi Science Support Cen- 
ter [32], using Pass 7 photon data. Our analysis uses 
Fermi Science Tools software version v9r27pl, released 
April 18, 2012. The first data selection are SOURCE-class 
photons from 200 MeV to 100 GeV in the region within 
5° radius of the origin of the Galactic coordinates. The 
maximum zenith angle is set to the Fermi-LAT recom- 
mended 100° to remove Earth limb effects, and the good 
time interval selection is set to the recommended selec- 
tion^] From this, we bin photons into 20 logarithmically- 
spaced energy bins in a 7° x 7° region of interest (ROI) 
square that fits within the initial selection circular region. 
For parameter and source fitting, we perform a binned 
likelihood analysis which generally follows the Extended 
Source Analysis Thread described at the Fermi Science 
Support Center [33]. The count maps for several energy 
bins are shown in the top row of Figs, [l] [2j & [3] The 
analysis procedure generates model maps from the model 
definitions of point and extended sources and varies pa- 
rameters to maximize the log-likelihood for the Poisson 
counts summed spatial and energy bins, defined as 



ln£ = ^ ki In — — In (fej), 



(3-1) 



where \±i is the model counts from a linear combination 
of all sources in the pixel z, and k{ is the observed counts 
in the pixel. Note that the Fermi Science tool gtlike 
reports the value for — ln£ — J2j lnfc Q- We quote the 
full — ln£ as computed from Eq. (3.1) in this paper. 



We generate the relevant 2FGL point sources that 
could contribute to the field of view using the user tool 
make2FGLxml .py [34]. The point sources included in our 
analysis are 17 point sources within the 7° x 7° ROI, 35 
sources beyond the 7° x 7° square region, two extended 
sources W28 and W30 associated with supernova rem- 
nants, the extended source 2-year Galactic diffuse map, 
and the diffuse isotropic component. The 17 sources in 
the ROI are varied in amplitude and spectrum, unless 
their point source test statistic (TS) significance is less 
than unity. In this case, the amplitude and spectrum 
of the source is fixed. The source identified with Sgr 
A*, 2FGL J1745.6-2858, was found to be better fit by 
a log-parabola than a broken power law, and therefore 
a log-parabola for its spectrum was chosen. This differs 
from the results of Refs. [22 , 35 which found a broken 
power law as a better fit for Sgr A*, though our obser- 
vation period contains a significantly larger time period 
than those studies. The log-parabola was also the spec- 
trum designation by the 2FGL catalog. The quantity 



1 (DATA_QUAL==l)&&(LAT_C0NFIG==l)&&(abs(R0CK_ANGLE)<52) 



TS is defined as twice the difference between the log- 
likelihood with (ln£i) and without (ln£o) the source, 
i.e. TS = 2(ln£i - ln£ ) (TS = 25 corresponds to a 
detection significance of ~ 5<r) [36] . Some point sources 
contribute only significantly below 200 MeV, therefore 
their significance drops in our > 200 MeV analysis. If the 
TS for the point source within the ROI is found to be be- 
low 25, the spectrum of the source is fixed. There is only 
one point source within 2° radius that has a TS that falls 
below 25 in our analyses, and that is 2FGL J1754. 1-2930, 
but our analyses and conclusions are not affected by fix- 
ing or varying this source. The 35 point sources and two 
extended sources, W28 and W30, beyond the ROI are left 
fixed to their 2FGL parameters (not varied in the likeli- 
hood analyses), but may contribute photons to the region 
due to the large point spread function of sources, partic- 
ularly at low energy, and are therefore included in our 
model generation. During the preparation of this report, 
two new point sources were identified near the GC in 
Ref. [19] : one previously included in the Fermi-LAT Col- 
laboration First Year Gamma-ray LAT Catalog (1FGL), 
1FGL J1744.0-2931c, and a new source named 'bkgA' by 
Ref. [19 . To test the effect of these two new sources 
on our analysis, we added them to the baseline model 
and the model with the best fit extended source model 
(7 = 1.2 and a log-parabola extended source spectrum). 
When added to the baseline model, 1FGL J1744.0-2931c 
was found with TS = 140.5, and bkgA was found with 
significance of TS = 158.7. When added to the best fit 
extended source model, the TS of the extended source 
was reduced by 20%. The extended source is still found 
at very high significance TS = 3371.9. Since they do not 
significantly affect the extended source results, we do not 
include these two new point sources in the other > 200 
MeV analysis runs. We include them both in all the runs 
with > 1 GeV cut. Both these analyses are described in 
detail below. 

We discuss the more complicated > 200 MeV energy 
cut analysis. Since the GC ROI is so crowded, and the 
sources' fluxes and spectra may have changed since the 
generation of the 2FGL point source parameter fits, we 
refit the point source flux amplitudes and spectra in the 
3.8 year data using the python-based binned likelihood 
Fermi-LAT Science Tools. In order to find a robust fit 
to the region, we start a baseline fit to the ROI using 
only the 2FGL point and known extended sources. The 
source parameters are allowed to vary sequentially from 
their distance to the center of the ROI (GC), initially 
fitting to only the amplitudes of sources within 2° radius, 
then to the full spectral model within 2° radius. Then, 
the remaining sources' amplitudes within the full 7° x 7° 
ROI are varied, and finally, all sources within the ROI 
and TS > 25 have their spectra varied. 

The analysis performed in the above paragraph is ini- 
tially done for only the known 2FGL point and extended 
Galactic diffuse amplitude and isotropic amplitude. This 
is our "baseline" model. 

We also perform a second analysis keeping only higher 
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1.76 - 2.40 GeV 2.40 - 3.28 GeV 3.28 - 4.47 GeV 4.47 - 6.10 GeV 




FIG. 2. Shown in the top row are photon counts in four energy bins that have significant evidence for an extended source with 
a spectrum, morphology, and rate consistent with a 100 GeV mass WIMP annihilating to 66-quarks in the 7° x 7° region about 
the GC. The panels show fits and residuals in the same manner as Fig.[l] but for higher energies. The maps have been filtered 
with a Gaussian of width a = 0.3°. The 17 point sources in the ROI are marked as circles in the top panels. 
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energy photons from 1 — 100 GeV in 8 logarithmically 
spaced bins on a slightly larger 10° x 10° region around 
the Galactic Center and with finer spatial binning of 0.05 
degree. As highlighted by Ref. [19], new point sources 
become significant in this energy band and therefore we 
check that our results are robust to a change in the spec- 
trum of photons analyzed. To enable direct comparison 
to the recent results from Ref. [19] , we only keep photons 
from Aug-4-2008 to Aug-4-2011. We vary sources within 
the inner 2 degrees and some other significant sources 
to converge to the baseline (high energy) model. The 
two new point sources found by Ref. [19] are included in 
the best-fitting models for this analysis, but not in the 
baseline model. Due to the lack of lower energy photons, 
convergence is more easily achieved as opposed the case 
where we include photons down to 200 MeV. 

To test the presence of any new extended source in the 
GC, we generate a number of extended source templates. 

• A profile with projected density index T = 0.7 [37] 
that is consistent with the stellar density profile of 
the nuclear stellar cluster. Note, however, that the 
bulk of the extended emission originates from out- 
side the region where Ref. [37] estimate the stellar 
density profile. 

• A set of seven p 2 templates (labeled "Density 2 " in 
the tables) with p chosen to be centrally-pe aked : 
six that are derived from a/3j profiles, Eq. (2.1) 
with a = 1, = 3 and 7 = 0.9, 1.0, 1.1, 1.2, 1.3, 1.4. 
The inner-profile slope 7 is the primary determi- 
nant of the signal morphology in the GC. However, 
in order to map our results on to the dark mat- 
ter annihilation cross section and particle mass pa- 
rameter space, we need to consider the full profile. 
The seve nth profile we adopt is an Einasto profile, 
Eq. (2.2), as an example of the prediction of dark 



matter only simulations. 

• To test for a dependence on the spatial morphol- 
ogy of the extended source, we also consider ax- 
isymmetric projected density profiles with axis ra- 
tio of 1:2 (labeled "Axisym" in the Table [Tv]) for 
the 1 - 100 GeV analysis with T = 0.7 and T = 1.4. 
We motivate the choice of T = 1.4 in ^V) 

Since the nature of the extended emission is uncertain, 
we adopt several spectral models for the extended emis- 
sion, including general log-parabola spectra, 



dN 
dE 



-{a+p \og(E/E b )) 



(3.2) 



with two parameters a and /3, and where is an arbi- 
trary fixed scale energy. We also test an extended source 
spectrum power-law with exponential cut-off, 



dN 
dE 



Nn 



E 



TABLE I. The best-fit TS, negative log likelihoods, and AC 
from the baseline for general models in the 200MeV - 100 
GeV analysis. 



Spatial Model 


Spectrum 


TS 


-InC 


Aln£ 


Baseline 






14(J(J7(J.Z 




Density r = 0.7 


LogPar 


1725.5 


139755.5 


314.7 


Density 2 7 = 0.9 


LogPar 


1212.8 


139740.0 


330.2 


Density 2 7 = 1.0 


LogPar 


1441.8 


139673.3 


396.9 


Density 2 7 = 1.1 


LogPar 


2060.5 


139651.8 


418.3 


Density 2 7 = 1.2 


LogPar 


4044.9 


139650.9 


419.2 


Density 2 7 = 1.3 


LogPar 


7614.2 


139686.8 


383.4 


Density 2 Einasto 


LogPar 


1301.3 


139695.7 


374.4 


Density 2 7 = 1.2 


PLCut 


3452.5 


139663.2 


407.0 


TABLE II. The best-fit TS, 


negative log likelihoods, and AC 



from the baseline, for specific dark matter channel models, 



using the a/3j profile (Eq. 2.1) with a 



3,7= 1.2. 



channel, m x 


TS 


-ln£ 


Aln£ 


bb, 10 GeV 


2385.7 


139913.6 


156.5 


bb, 30 GeV 


3460.3 


139658.3 


411.8 


bb, 100 GeV 


1303.1 


139881.1 


189.0 


bb, 300 GeV 


229.4 


140056.6 


13.5 


bb, 1 TeV 


25.5 


140108.2 


-38.0 


bb, 2.5 TeV 


7.6 


140114.2 


-44.0 


r + r~, 10 GeV 


1628.7 


139787.7 


282.5 


t + t~, 30 GeV 


232.7 


140055.9 


14.2 


t + t~, 100 GeV 


4.10 


140113.4 


-43.3 



~{E/E C ) 



(3.3) 



with power law 7, cut-off energy E c and arbitrary fixed 
scale energy Eq. 

For the dark matter halo models, we also include spec- 
tra of photons from dark matter particle annihilation 
into bb quarks and r + r~ leptons for dark matter particle 
masses of 10, 30, 100, 300, 1000 and 2500 GeV, gener- 
ated with pythia 6.4 as described in Ref. [38]. Nearly ev- 
ery combination of morphology and spectrum was walked 
through the iterative parameter relaxation procedure de- 
scribed above. 



IV. RESULTS 

The iterative fitting procedure described in the previ- 
ous section revealed significant detections of an extended 
source in the GC. The model fits found numerically con- 
vergent fits for several spatially extended sources with 
a number of spectral types. Importantly, the extended 
source has a strong degeneracy with the several point 
sources nearest the GC. The four point sources nearest 
the GC, Sgr A* (2FGL J1745.6-2858), 2FGL J1746.6- 
2851c, 2FGL J1747.3-2825c and 2FGL J1748.6-2913 re- 
duce the amplitude of their emission from the baseline 
model to the 7 = 1.2 Density 2 Log-Parabola spectrum 
model by factors of 3.1, 1.21, 1.9, 2.0, respectively. This 
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FIG. 3. Shown in the top row are the photon counts in four energy bins in the 7° x 7° region about the GC that have significant 
evidence for an extended source with a central morphology consistent a projected density-squared map with a central density 
profile index 7 = 1.3. This could be consistent with a concentrated population of unresolved points sources as discussed in 
j ]V[ The extended source is best-fit with a Log-Parabola spectrum. The panels show fits and residuals in the same manner as 
Fig.[l] The maps have been filtered with a Gaussian of width a — 0.3°. The 17 point sources in the ROI are marked as circles 
in the top panels. 
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TABLE III. The best-fit total flux and 68% error of the GC 
extended source models for the 200 Me V- 100 GeV analysis. 
LP is log-parabola spectrum, and PLcut is power-law spec- 
trum with an exponential cut off. 



model 



Density r = 0.7 LP 



Density 2 
Density 2 
Density 2 
Density 2 
Density 2 
Density 2 
Density 2 
7 1.2 
7 1.2 
7 1.2 



7 = 0.9 LP 
7 1.0 LP 
7 1.1 LP 
7 1.2 LP 
7 = 1.3 LP 
Einasto LP 
7 = 1.2 PLcut 
66, 30 GeV 
66, 100 GeV 
t + t~ 10 GeV 



flux and error [ph cm s ] 



(1.31 ±0 
(2.31 ±0 
(5.29 ±0 
(3.36 ± 
(2.69 ±0 
(2.01 ±0 
(4.21 ±0 
(2.97 ±0 
(1.77 ±0 
(4.90 ± 
(5.13 ±0 



.06) x 10" 5 
,.06) x 10" 6 
).40) x 10" 6 
.23) x 10" 6 
.17) x 10" 6 
.11) x 10" 6 
,.32) x 10" 6 
1.22) x 10" 6 
).06) x 10" 6 
.23) x 10" 7 
.20) x 10" 7 



TABLE IV. The best-fit negative log likelihoods, AC from 
the baseline model and fluxes with 68% errors for the general 
models in the 1 - 100 GeV analysis. The baseline model 
for this analysis has ln£ = —176478.6. LP is log-parabola 
spectrum, PLcut is power-law spectrum with an exponential 
cut off and PL is power-law spectrum without an exponential 
cut off. 



Spatial Model 


Spectrum 


Aln£ 


flux [10 -7 ph cm" 2 s" 1 ] 


Density 2 7 = 1.0 


LogPar 


189.5 


1.57 ±0.08 


Density 2 7 = 1.2 


LogPar 


206.2 


1.51 ±0.09 


Density 2 7 = 1.2 


PLcut 


205.4 


1.49 ±0.09 


Density 2 7 = 1.2 


PL 


126.1 


1.22 ±0.08 


Density 2 Einasto 


LogPar 


189.2 


1.45 ±0.09 


Axisym r = 1.4 


LogPar 


202.1 


2.00 ±0.12 


Axisym r = 0.7 


LogPar 


165.5 


1.87 ±0.15 



indicates that the central point source fluxes are in- 
creased by the baseline 2FGL model in order to try to fit 
the presence of the extended emission. The central point 
sources' spectra change significantly as well. Since point 
sources are over-subtracting the extended source, it leads 
to the appearance of "holes" in the emission residuals of 
the extended source, as seen in the second row of Figs.[l} 
[3J Note that there is an oversubtraction near the position 
of (6, i) = (0, -2°) (or 0, 358°) that is due to a feature at 
that position in the Galactic diffuse model. 

The extended sources were found to be detected at 
high TS significance for several spectral models. The TS 
magnitude and best-fit log likelihood values to several 
general morphological models and spectra are shown in 
Table H Note that the TS differs from 2 A In C because 
all other model components are also changing, and the 
presence of the source is not the only change relative to 
the baseline model. This also implies that while the TS 
has a definite statistical interpretation, its significance 
once it becomes large is questionable. Results to fits with 
several dark matter particle mass cases and annihilation 



channels to 66 quarks and r + r~ channels are given in 
Table [nj There is a good fit for the extended source model 
with dark matter particles with masses from 10 GeV to 
1 TeV annihilating into bb quarks, and particle masses of 
10 GeV to 30 GeV annihilating into r + r~ leptons. As 
shown in Table [il] the upper limit of the particle mass 
to give a significant detection, TS > 25, is between 1 
and 2.5 TeV in the case of annihilation to bb quarks, 
and between 30 and 100 GeV in the case of annihilation 
into r + r~ leptons. This is a significant finding and hints 
at the possibility of an underlying signal due to dark 
matter even if the bulk of the extended emission is due 
to astrophysical sources. 

The best-fit model for an extended source in the GC is 
a projected density-squared source with 7 = 1.2 for the 
inner density profile and a general log-parabola spectrum. 
The spectrum is best-fit by a log-parabola with 7V = 
(3.17±0.33) x 10- 3 ph cm" 2 s" 1 sr" 1 , a = 0.488±0.062 
and p = 0.325 ± 0.011, with fixed E b = 100 MeV. The 
spectrum of the extended source is also consistent with 
a power-law with exponential cutoff with Nq = (6.62 ± 
0.74) x 10~ 3 ph cm" 2 s" 1 sr" 1 , a = 1.48 ± 0.05 and 
E c = 2.46 ±0.2 GeV, with E = 100 MeV. Other central 
profile index values for 7 as well as the Einasto profile 
gave good fits and were detected at high TS significance. 
The case of 7 = 1.4 was found to be strongly degenerate 
between the extended source spectrum and amplitude 
and that of Sgr A* (2FGL J1745.6-2858), since the flux 
of the extended source was largely within the PSF of the 
Fermi-LAT spatial resolution. 

Since the dark matter particle mass is a prior for the 
GC gamma-ray analysis, we treat it as a systematic un- 
certainty. The best-fit dark matter annihilation models 
we tested are the case of dark matter particle masses, m x 
of 30 GeV and 100 GeV annihilating to bb and a a/3j dar k 
matter profile with a = l,/3 = 3,7 = 1.2 (cf, Eq. [2T| . 
Using the model of m x = 30 GeV annihilating into bb, 
shown in Fig[l]are: the data counts map (first row), base- 
line residuals (second row), GC extended source model 
best fit (third row), residuals when not including the ex- 
tended source in the best fit model (fourth row), and 
total model residuals (bottom row), for four significant 
energy bins. Because the dark matter particle mass is a 
systematic uncertainty, though the m x = 100 GeV case 
is not as good of a fit as measured by the A In £, this case 
is significant and still detected at TS significance. Tak- 
ing m x = 100 GeV annihilating into 66, the data counts 
map (first row), baseline model residuals (second row), 
GC extended source model best fit (third row), residu- 
als when not including the extended source in the best 
fit model (fourth row), and full model residuals (bottom 
row), for four significant energy bins, are shown in Fig[5J 

The 7 = 1.3 density profile model data counts map 
(first row), baseline residuals (second row), GC extended 
source model best fit (third row), residuals when not in- 
cluding the extended source in the best fit model (fourth 
row), and total model residuals (bottom row), for four 
significant energy bins, are shown in Fig [3] The best- 
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FIG. 4. Shown the parameters of particle dark matter mass 
m x and cross section (av) for annihilation to bb quarks consis- 
tent with the extended gamma-ray source in the GC at 68% 
CL (dark pink) for a d ark matter density profile with central 
slope 7 = 1.2 (cf, Eq. 2.1), our best fit spatial model. The 
red line is the case of pq — 0.3 GeV cm -3 . The diagonally 
hatched region is approximately where the TS significance 
drops below ~ 5cr. The light pink region shows the extension 
of the consistency region for 7 = 1.3, with vertically hatched 
region corresponding to approximately where the TS signifi- 
cance drops below « 5a. The region above the solid line in- 
dicates the parameters excluded at 95% CL by stacked dwarf 
analyses ^0]. The region above the dashed line indicates the 
parameters excluded at 95% CL by HESS observations of the 
GC [12] . We have assumed here that all of the extended emis- 
sion is due to dark matter annihilation. If only part of it is 
due to dark matter, then the required cross section should be 
lower. 




100 



m x [GeV] 



FIG. 5. Shown are the parameters of particle dark matter 
mass m x and cross section (av) for annihilation to r + r~ lep- 
tons consistent with the extended gamma-ray source in the 
GC at 68% CL for a central density profile of 7 = 1.2 (the 
best-fit model, in dark pink) and 7 = 1.3 (light pink). The red 
line is for the case of p© = 0.3 GeV cm -3 . The diagonally 
and vertically hatched regions are approximately where the 
TS significance drops below « 5cr for the 7 = 1.2 and 7 = 1.3 
cases, respectively. The region above the solid line indicates 
the parameters excluded at 95% CL by stacked dwarf analy- 
ses [T0l. 



exclusion bounds from the stacked MW satellite analysis 
(see 



DISCUSSION 



fit spectrum for this model is a log-parabola with Nq = 
(2.33 ± 0.39) x 10- 3 ph cm" 2 s" 1 sr~\ a = 0.47 ± 0.11 
and p = 0.328 ± 0.019, with fixed E b = 100 MeV. 

The results from the 1 — 100 GeV analysis lends further 
support for the results described above. The differences 
in Aln£ values for Einasto and 7 = 1 models vs the 
best-fit 7 = 1.2 model are still visible in the 1 — 100 GeV 
analysis. Thus, the higher energy analysis also seems 
to prefer a projected density-squared map with 7 = 1.2. 
We also allowed the projected density maps to be axisym- 
metric with axis ratio 1:2. There was no significant dif- 
ference between the annihilation model (Density 2 ) with 
7 = 1.2 and the axisymmetric projected density model 
with T = 1.4 (both with LogParabola spectrum). Thus 
this preliminary analysis of the 1 — 100 GeV data indi- 
cates that the data is unable to pick out a morphology 
for the extended emission. 

The total flux in the 1 — 100 GeV range for the Density 2 
7 = 1.2 LogParabola spectrum model is about a factor 
of 2.7 smaller. Our main analysis for the limits on dark 
matter particle mass and cross section was performed 
with the > 200 MeV cut. Thus, if instead, we were to use 
the > 1 GeV cut, the required annihilation cross section 
will be lower and this will decrease the tension with the 



There is definitive evidence from our analysis that 
there exists a new source in the Galactic Center that 
is not associated with any sources within the 2FGL or 
Fermi-LAT Galactic diffuse maps. Below we discuss some 
interpretations of the results. 

A. Dark Matter Annihilation 

Significantly, we find a good fit when using gamma- 
ray spectra arising from dark matter annihilation. The 
fits are consistent in morphology, spectrum, and, as we 
show, in annihilation rate expected in thermal dark mat- 
ter production models. The fits are consistent with a 
wide range of particle mass annihilation spectra, from 
from 10 GeV to 1 TeV dark matter particles annihilating 
into bb quarks , as well as from 10 GeV to 30 GeV par- 
ticle mass annihilating into r + r~ leptons. We have not 
performed an exhaustive search of the parameter space 
of relative annihilation channels, particle mass, and halo 
morphology. 

We explore the parameter space consistent with the 
dark matter interpretation by varying the primary un- 
certainties in the signal: the scale density p s , which 
is set by the local dark matter density p = 0.3 ± 
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0.1 GeV cm -3 [30[|3T], and inner halo profile 7. Since 
the signal is proportional to the density- squared integral, 
which is normalized by p si the error propagation leaves 
the uncertainty in J in any direction to be double that in 
p s , assuming Gaussian errors as an approximation. The 
resulting uncertainty in the annihilation rate is commen- 
surate to that in J, and therefore a wide range of annihi- 
lation rates are consistent with the signal. The range of 
particle annihilation rates (av) into bb quarks and dark 
matter particle masses m x consistent with the signal are 
shown in Fig. [4] The range of annihilation rates (av) 
into r + r~ consistent with the signal are shown in Fig. [5] 
Note that the solid bands in Figs. |4][5] are the ranges 
where annihilation into these channels is consistent with 
the extended emission at greater than 5a. The Aln£ 
values however prefer m x ~ 30 GeV for the bb channel 
and m x ~ 10 GeV, similar to that found in Ref. [16], 
though the central best-fit region of our fits prefer higher 
values of annihilation cross-section, largely due to the 
higher flux attributed to the extended source relative to 
non-Sgr A* point sources, which were fixed in that work. 

Due to the significant uncertainty in the amplitude of 
the dark matter signal source J, there is a large range of 
parameters consistent with the source in the GC. Parts 
of the parameter space have been excluded at 95% CL by 
stacked dwarf analyses [H[T0], and, in the case of the bb 
channel, are bordering on that excluded by observations 
of the GC by HESS [12]. In this regard, we note that the 
1 — 100 GeV energy cut analysis prefers a lower annihila- 
tion cross section and hence the tension with the results 
from the stacked dwarf analysis could be weaker. How- 
ever, there are parts of the parameter space that are still 
consistent with all other constraints, and, most signifi- 
cantly, have an annihilation cross section in the vicinity 
of thermal relic value (av) ~ 3 x 10 -26 cm 3 s _1 at a typ- 
ical WIMP mass of ~ 100 GeV. Prior work has found 
that the GC extended source is consistent with a nar- 
rower range of parameters, with annihilation in the bb 
channel at a mass-scale of 30 GeV or into r + r~ with a 
mass scale of 10 GeV, with both cases having a narrow 
(av) w 1.0 x 10~ 26 cm 3 s" 1 [16]. 



B. Pulsars and other point sources 

The apparent extended morphology can be a simple 
superposition of several point sources close to the GC in 
projection. The importance of proper point source sub- 
traction has been emphasized in Ref. [22]. As described 
in j |III| two new point sources were found in Ref. [19] 
in the GC. A large number of unresolved point sources 
can reproduce the features of an extended source if their 
3D density profile falls off steeply enough. Here we con- 
sider whether millisecond pulsars (MSPs) have the right 
properties to explain the extended emission. 

The spectrum of gamma-ray emission from MSPs in 
the GC would be consistent with that expected from 
stellar globular clusters such as Omega Cen and 47 Tuc 



[T7J [39] and this in turn is consistent with the log- 
parabola fits we find here. The values of Aln£ favor 
a compact object population that follows a 3D number 
density profile n(r) with a log-slope — <iln(n)/<iln(r) = 
27 = 2.4. To see if this would be consistent with an 
unresolved MSP population in the Galactic Center, we 
look at resolved low mass X-ray binary (LMXB) popula- 
tions, which should have a similar spatial profile as that 
of MSPs. Indeed, LMXBs and MSPs are thought to be 
different phases of the same binary system. Observations 
targeting LMXBs in M31 show a sharp rise in the sur- 
face density within about an arcminute [40] . Outside this 
region, the LMXBs track the K-band luminosity. The in- 
ner "excess" is consistent with a population created by 
stellar encounters in the extremely high density environ- 
ment in the central regions of the Galactic bulge [41]. 
The physical scale corresponding to 1 arcminute is about 
200 pc, which at the distance of the MW center would 
be approximately 1.5 degrees. This is exactly the region 
from which the bulk of the extended source emission em- 
anates. We estimate a power-law index of —1.5 ± 0.2 for 
the projected M31 LMXB distribution [40] between 10 
and 100 arcmin. The projected distribution correspond- 
ing to our best-fit LogParabola spectrum density-squared 
model (which has 7 = 1.2 is R~ 1A (where R is the pro- 
jected radius), consistent with the surface density profile 
of the inner M31 LMXB population. 

The LMXB population in the center of the Milky Way 
is less well determined. A study using INTEGRAL [42] 
found too few LMXBs in the inner 1 degree radius to 
robustly infer a profile but there was slight evidence of 
steepening compared to the stellar profile in the tran- 
sient LMXBs that may be consistent with the dynamical 
formation scenario [41]. Thus, both the Milky Way and 
the M31 LMXB population comparisons lend support to 
our proposal (to different degrees) that the spatial distri- 
bution of gamma-ray bright stellar remnants in the GC 
could be steeper than 1/r 2 . 

The flux in the extended source for the density squared 
model with 7 = 1.3 or a projected ma p of 1/r 



2.6 



is 



(2.01 ±0.11) x 10~ 6 ph cm" 2 s" 1 (cf, Table III). It would 



interesting to estimate how many MSPs would be re- 
quired to account for all of this emission. We choose 
47 Tuc as a reference. The measured flux between 0.1 
and 10 GeV is 2.8(±0.6) x 10- 8 ph cm" 2 s" 1 [43 . Us- 
ing the best-fit model and assuming a population of 30- 
60 MSPs in 47 Tuc gives us a typical GC MSP flux of 



1 - 2 x 10" y ph cm" 



in the 0.2 to 100 GeV band. 



Thus we see that ~ 1000 MSPs are required in the in- 
ner bulge region to explain all of this extended emission. 
This conclusion is unchanged if we adopt the measured 
flux for MSPs in our Galaxy unassociated with globular 
clusters [44] . 

To get a sense for how reasonable this hypothesis is, 
we compare the required number of MSPs to the stel- 
lar content with the inner bulge region. Within the 
dense molecular clouds of the central few hundred par- 
sees ( "central molecular zone" ) [45] [46] , there is a com- 



11 



pact region named the "nuclear bulge" - projected radius 
smaller than about 2 degrees - that seems to be distinct 
from the old Galactic bulge population [46j [47]. Esti- 
mates of the stellar content based on the near infra-red 
luminosity suggest a total stellar mass of ~ 1O 9 M [47] 
and most this mass is within the inner 1 degree. This is 
^ 1000 times more than the stellar mass in 47 Tuc glob- 
ular cluster and the required number of MSPs is about 
~ 20 times more than that in 47 Tuc. This is plausi- 
ble despite the large velocity dispersion in the Galactic 
Center given the higher stellar densities in the Galactic 
nuclear bulge. Putting these observations together with 
the suggestive M31 LMXB steep inner density profile al- 
lows us to make the reasonable argument that the bulk of 
the extended emission in the ~ GeV energy range could 
be is associated with the stellar content of the nuclear 
bulge. 

A stellar projected density profile of T = 0.7 is also 
consistent with the emission, though less preferred (cf. 
Table [I]). The spectrum of the emission from these source 
is consistent with either log-parabola or a power-law with 
an exponential cutoff, though a log-parabola is favored 
(Aln£ = 9.3 between these two models). A good fit 
is also achieved by a power-law with exponential cutoff 
spectrum as in Eq. (3.3), which is expected from MSPs 



as those known to exist in globular clusters [17], though 
the scale of the exponential cutoff is slightly higher (E c = 
2.46 ± 0.2 GeV) than that observed for globular clusters 
(E c w (1.0 to 1.8 GeV) ± 1 GeV), but not significantly 
so given errors on the globular cluster spectra. 



C. High-energy Cosmic Rays Interacting with Gas 

The GC source may also be consistent with gamma-ray 
emission from cosmic rays interacting with gas within the 
inner 3 pc to 300 pc of the region near Sgr A* [T6 | fT8| fT9] . 
Two mechanisms have been proposed: (1) from cosmic 
ray protons on gas within the inner ~3 pc leading to 
hadronic p — p collision gamma-rays [18] and (2) cosmic- 
ray electrons producing bremsstrahlung gamma-rays on 
molecular gas [19 . In the case of hadronic emission, the 
flux has been found to be extended but within the PSF 
of the Fermi-LAT [18 . Therefore, though this could be 
a contribution to the emission in the GC region, it does 
not account for the significant evidence for an extended 
source. 

In the case of cosmic-ray electrons producing gamma- 
rays via bremsstrahlung on the molecular gas, there can 
be a significant spatial extent to the emission. Ref. [19] 
finds that the source electron population is consistent 
with radio observations of synchrotron emission from the 
high-energy population of electrons, as well as the mor- 
phology of the Fel 6.4 keV X-ray emission. In addition, 
they find that using the radio emission morphology, trac- 
ing the synchrotron emission from the cosmic ray elec- 
trons, improves the fit to the observed extended gamma- 
ray emission by 2Aln£ = 113, and the observations are 



consistent with the model's energy spectrum from 1 GeV 
to 100 GeV. Our 1 — 100 GeV analysis mirrors that of 
Ref. [19 in pixel resolution and ROI and the time pe- 
riod was chosen to be the same for the purpose of com- 
parison. Thus it is worth noting that the improvement 
we obtain for the Density 2 7 = 1.2 LogParabola model 
is 2Aln£ = 412, significantly better than that obtained 
using the 20 cm radio emission template. Our 7 = 1.2 PL 
model with only the galactic diffuse, isotropic, extended 
source and Sgr A parameters (8 in all) varied, and not 
including the two new sources in Ref. [19 , should be a 
better comparison to the radio emission template model. 
For this model, we obtained 2 A In C = 252 - a poorer fit 
compared to 7 = 1.2 model with the log-parabola or the 
PLcut (power-law with-exponential cut off) spectra, but 
a better fit than the radio emission template model. This 
clearly deserves further study but is beyond the scope of 
the present work. 



VI. CONCLUSIONS 

Our analysis has revealed a source in the Galactic Cen- 
ter at high significance that is consistent with extended 
emission. The most intriguing aspect of this source is 
its consistency in morphology, spectrum and flux with 
that expected from canonical thermal weak-scale parti- 
cle dark matter in a centrally peaked halo density pro- 
file. The best-fitting dark matter models have particle 
masses around 30 GeV annihilating to bb with a halo 
density profile that is somewhat steeper than the cold 
dark matter simulation predictions, consistent with the 
results of Refs. [I5j [16]. The source is also consistent 
with extended emission from a stellar remnant popula- 
tion or from bremsstrahlung of cosmic rays (produced 
around Sgr- A*) on molecular gas. Because the spectrum 
and rate of an astrophysical source interpretation is less 
well specified, a broader range of spectra and fluxes can 
be accommodated. The log-parabola and power-law with 
exponential cutoff spectra expected in these interpreta- 
tions are consistent with the observations. 

Occam's razor would dictate a conservative interpreta- 
tion of these results that strongly prefers the astrophys- 
ical explanations of the source signal. The bulk of the 
emission seen here is likely to be another piece in the 
puzzle of the violent processes involved in the crowded 
region near Sgr- A*, associated with cosmic ray interac- 
tions with molecular gas in the central 300 pc [19 , and 
from a centrally-concentrated MSP population [T7j . 

However, since the Galactic Center is also the region 
with the highest expected luminosity in gamma-rays due 
to dark matter annihilation, the three-fold consistency 
of morphology, spectrum and rate with that which is 
expected from canonical weak-scale thermal dark mat- 
ter should not be dismissed. Our results confirm that 
of Refs. [I5j [16] in finding significant evidence of an ex- 
tended source in the GC, but we find that a broader set of 
source spectra, dark matter particle masses and annihila- 
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tion rates are compatible with the data. This is primarily 
because the spatial response of Fermi-LAT changes with 
energy, and the complex crowded region requires a si- 
multaneous fit of point sources, diffuse emission as well 
as any new extended source morphology and spectrum. 
This results in a much broader consistent model space 
than single-region fixed- astrophysical-source spectral fits. 
The dark matter interpretation of the gamma-ray signal 
can be complicated by the existence of the other poten- 
tial extended sources in the GC, and the flux from dark 
matter may be lower than our single-extended-source fits 
provide. This would prefer lower annihilation cross sec- 
tions than our single-component models find. 

Further measurements toward dwarf galaxies, the 
Milky Way halo, or simultaneous analyses of multiple- 
regions could reach significantly into the parameter space 
consistent with the dark matter interpretation. It would 
take indirect detections towards multiple sources with 
equivalent spectra, particle dark matter mass, and an- 
nihilation rates to affirm a beyond the standard model 
interpretation of the source in the Galactic Center. 
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